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ABSTRACT 

If a galaxy cluster is effectively in dynamical equilibrium then all galaxy populations 
within the cluster must have distributions in velocity and position that individually 
reflect the same underlying mass distribution, although the derived virial masses can 
be quite different. Specifically, within the CNOC cluster sample the virial radius of 
the red galaxy population is, on the average, a factor of 2.05 ± 0.34 smaller than that 
of the blue population. The red galaxies also have a smaller RMS velocity dispersion, 
a factor of 1.31 ± 0.13 within our sample. Consequently, the virial mass calculated 
from the blue galaxies is 3.5 ± 1.3 times larger than from the red galaxies. However, 
applying the Jeans equation of stellar-hydrodynamical equilibrium to the red and 
blue subsamples separately give statistically identical cluster mass profiles. This is 
strong evidence that these clusters are effectively equilibrium systems, and therefore 
empirically demonstrates that the masses in the virialized region are reliably estimated 
using dynamical techniques. 

Subject headings: galaxies: clusters, cosmology: large-scale structure of universe 

1. Introduction 

The primary goal of the CNOC (Canadian Network for Observational Cosmology) cluster 
redshift survey is to obtain a value of density parameter, Qq, for those components of the mass 
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field (possibly all) that participate in gravitational clustering (Carlberg et al. 1996, hereafter 
CNOCi). The product of the field luminosity density, j, with its mass-to-light ratio, M/L, 
estimates the mean mass density of the universe, po ( Port 1958| ) , which in ratio to the critical 
density, p c is equal to Qq ( pott fc Turner 1976j ), when calculated in current epoch co-ordinates. 
The crucial quantity is an estimator of the field mass-todight ratio, for which we use the cluster 
virial mass-todight ratio, M v /L. The main technical concern of the CNOC survey is that there are 
various possibilities that M v is a biased estimator of the cluster mass. Besides statistical biases, 
dissipative processes can cause the cluster galaxies to become more centrally concentrated than 
the cluster mass. In addition, when estimating the field M/L from the cluster results, allowance 
must be made for the differences between the cluster and field galaxy populations. 

To address these issues a sample of high X-ray luminosity galaxy clusters were selected from 
the EMSS catalogue ( [Henry et al. 1992j , |G ioia &: Luppino 1994] ) . This produces a relatively 



homogenous set of 16 clusters that likely contain a substantial virialized component. The sample, 
at a mean redshift of about 73, gives a substantial column in redshift space which we use to 
calculate the field luminosity density in precisely the same measurement system as the cluster 
light. The galaxies were selected with no knowledge of the colour or whether they were cluster 
members, although great care is taken to control the selection effects that are unavoidable to 
maintain efficiency ( Yee, Ellingson Carlberg 1996| ). 



From these data we calculated the line-of-sight velocity dispersions, cxi, and virial radii, t V: 
of the clusters to derive the virial masses ( Binney fc Tremaine 1987 ), GM V = 2>a\r v . The velocity 



dispersions are calculated with an iterative technique that statistically removes interlopers by 
taking advantage of the relatively large foreground and background sample to estimate interloper 
densities ( pNOCi ). This technique leads to velocity dispersions that are on the average 7% less 



than those found from precisely the same data in the redshift range, and about 13% less with a 
somewhat less restrictive redshift cut, even when evaluated with the iterated bi- weight estimator 
(Beers, Flynn & Gebhardt 1990). The virial radius is calculated from the data in the redshift 



range of the cluster with a "ring-wise" potential estimator, which allows correction for the strip 
sampling and helps reduce the noise (details are given in CNOCi| ). All these global quantities 



are derived from the individual galaxy positions, velocities, colors and luminosities, for which 
objective errors are calculated using the Jackknife method ( |Efron &: Tibshirani 1986] ). The errors 
are sufficiently small that the assumption of a symmetric distribution is appropriate and is verified 
by comparing to Bootstrap error estimates. 



The average cluster virial mass-to-light (r band, k-corrected) ratio is (M v /L 



k\ 



289 ± 50/iMq/Lq (CNOCi). In the field, over the same redshift range, with the same sampling 
and corrections, the closure mass-to-light ratio is p c /j = 1136 ± I38/1M0/L0 (in co-moving units, 
for qo = 0.1). There are two additional corrections applied to (M v /L^). The mean luminosity of 
galaxies brighter than of —18.5 mag is 0.11 ± 0.05 lower in the cluster than the field and the 
analysis of the mass profile (Carlberg, Yee & Ellingson 1997, hereafter CNOCii) showed that M v 



needed to be multiplied by 0.82 ± 0.14, which we attribute to neglecting the surface term in the 
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virial mass estimator. Making these corrections leads to a redshift adjusted = 0.19 db 0.06. 

Our conclusion that the stellar-hydrodynamical (SHD) mass profile of the clusters averaged 
together closely follows the number density profile of the galaxies (CNOCii) rests on the 



assumption that the clusters are effectively in equilibrium such that the Jeans equation can be 
used to derive the mass profile from the radial and velocity distribution of any tracer population. 
In this paper we demonstrate that two radically different tracer populations are consistent with 
the same underlying mass profile, which we take as strong support of the equilibrium assumption. 
In the following section we describe the division of the sample into two completely disjoint blue 
and red subsamples, whose surface density and projected velocity dispersion profiles are measured 
and demonstrate the gross differences in the calculated virial masses. In Section 3 we derive the 
stellar-hydrodynamical potential generating mass profiles in which the galaxies are orbiting. The 
implications of the statistical equality of the two profiles are discussed in the final section. All 
calculations in this paper assume Hq = 100 kms -1 Mpc -1 and qo = 0.1, although the results are 
either independent or not sensitive to these choices. 



2. The Blue and Red Subsamples 

The CNOC cluster redshift survey catalogues QYee, Ellin gson fc Carlberg 1996] ) contain ~2600 
galaxies having Gunn r magnitudes, g — r colors, and redshifts. The sample is split at the color 
(g - r ) z = (g - r)[1.2 + 2.33(z - 0.3)]" 1 of 0.7 QCNOCj ) . This is approximately the color of an 



Sab galaxy, which roughly divides galaxies into spheroid dominated, "old" stellar systems and disk 
dominated, "young" stellar systems. Fourteen clusters are averaged together after removing two 
"binary" clusters. 

To a limit of M r fc = -18.5 mag, about 70% of the cluster galaxies fall into the red subsample. 
The methods of analysis are precisely the same as used previously for the full sample flCNOCii ). 



The velocities are normalized to the o\ of each cluster calculated from galaxies of all colors so that 
the two subsamples have a common reference. The RMS velocity dispersions of the red and blue 
galaxies are 0.97 db 0.04 and 1.27 db 0.11, respectively, in units where the velocities are normalized 
to the velocity dispersions of all the galaxies in each cluster. That is, the velocity dispersion of 
blue cluster galaxies is about 30% higher than for red galaxies ( [Rood et al. 1972] , Kent & Gunn 



1982] [Stein 1996Q . The projected radii are normalized to r2oo = V3&i/(WH(z)), the radius at 



which the mean interior overdensity in a cluster is 200p c (2:). The quantity H(z) is the Hubble 
constant at the redshift of the cluster. The virial radii, r v , are 2.58 db 0.27 and 5.28 db 0.73 for red 
and blue galaxies, respectively, in r2oo units. In ratio to the sample as a whole, the average virial 
masses are 0.66 ± 0.08 and 2.31 ± 0.42, respectively, for red and blue subsamples. The ratio of 
these two virial masses is 3.5 ± 1.3. As is well known, the red galaxies in clusters are much more 



centrally concentrated than the blue galaxies (e.g. Dressier 1980 ). Our data usefully illustrate 
that the galaxy color selection can make a very large difference to the estimated mass. In general, 
blue or emission line selected galaxies will always give higher mass estimates than absorption line 
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galaxies. 

The projected surface number density profiles, Ejv(-K) of the two subsamples are displayed in 
Figure 0. The profil es are fitted to the projection of the density function ( [Hernquist 1990| ), 

u(r) = . - (1) 
w r(r + a) 3 

by adjusting A and a to minimize the residuals. The fits give scale radii a of 0.56 ± 0.10 and 
1.82 db 0.27 for red and blue populations, respectively. For the sample as a whole a = 0.66 ± 0.09. 
In both cases the x 2 P er degree of freedom is about 0.8, with the errors estimated as the average 
of the upper and lower la confidence range from the bootstrap analysis. 

The observed projected velocity dispersion, a p (R), can be used to infer the real space radial 
velocity dispersion, oy(r), via, 

f°° R 2 r 

a p 2 (R)Z N (R) = / v{r)a 2 .(l - 0—)—===dr, (2) 

where (3 = 1— a 2 , /a 2 , is the velocity anisotropy parameter. The parameters of 

b + r 

are adjusted until the x 2 with the observed a p {R) are minimized. The data and the fits for 
two values of (3 are displayed in Figure ||. Again the errors are from the symmetrized bootstrap 
confidence range. 



The Mass Distribution 



The mass distribution in which the galaxies are in orbital equilibrium can be inferred with 
no assumption about the relative distribution of mass and galaxies from the projected velocity 
dispersion profile, a p (R), and the projected galaxy density distribution, Sjv(-R)- Remarkably, we 



found (CNOCii) that the integrated galaxy number density profile, L(r), of our full sample is 
statistically indistinguishable from the integrated mass profile, M(r). 



Mass profiles are derived from the fitted u(r) and a r 
Tremaine 1987), 

, . air din a 2 dlnv 

M SH D(r) = -^r -77-^ + 37— 
G a In r a In r 



using the Jeans equation (Binney & 



(4) 



We will call this mass the stellar hydrodynamical mass. The resulting MsHD{ r ) from the red 
and blue galaxies are shown in Figure |3| for /3 = and 1 /2- The quantity plotted is the virial 
mass-to-light bias, 

M SH D(r) L 



bMv{r) 



L(r) M v 



(5) 
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where L is an arbitrary normalization of the total light that cancels between numerator and 
denominator. The quantity M v is the normalization of the virial mass of the dimensionless sample 
as a whole, M v = 3di 2 f v , where 6\ = 1.05 db 0.04 and r v = 3.35 ± 0.32. If the virial mass calculated 
from the sample as a whole gave the correct M/L at all radii then bMv( r ) = 1 everywhere. If 
the tracer population follows the true mass profile but with some scale error, then buv will be a 
constant other than unity. Figure ||| shows that beyond about 0.3r2oo the mass profiles deduced 
from the two subsamples are identical within their errors (about 15 and 30% for red and blue 
mass profiles, respectively). Figure |I] shows the value of buv evaluated at r2oo> demonstrating that 
the blue and red subsamples give statistically identical values. Furthermore the value is always 
less than unity, as we found for the sample as a whole ( |CNOCi]| ) which has the implication that 
although the galaxy numbers fairly accurately trace the cluster mass profile the virial mass to light 
ratio is always an overestimate of the mean M/L inside r2oo- There is no statistically significant 
radial gradient of bM v ( r ) for r ~ r 2 oo- 

There are several important conclusions to be drawn from Figure [3|. First, it provides strong 
evidence that the blue and red galaxies are sufficiently in equilibrium with the cluster mass 
distribution that these dynamical methods work to recover the true mass profile. Second, the 
fitted scale radius of the mass profile (taken to be equal to that of the full sample) a = 0.66 ± 0.09, 
is slightly more extended than the red galaxies ( Tyson fc Fischer 1995| ) a = 0.56 ± 0.10, but is 
substantially more compact than the blue galaxy distribution, a = 1.82 db 0.27, all measured in r2oo 
units. The differences in scale radii imply that neither the blue nor the red galaxies are distributed 
like the mass, so the virial masses calculated from these subsamples are not correct. 



4. Conclusions 

The densities and velocity dispersions of the red and blue galaxy subsamples, although very 
different from each other, imply a statistically identical mass profile for the cluster potential. This 
is necessary if both populations are in equilibrium with the potential, validating our use of the 
Jeans equation. This mass profile is statistically identical to the number density distribution of 
all the galaxies of our full sample. The limiting factor in the precise numerical agreement is that 
there are relatively few blue galaxies in clusters, meaning that their density and velocity dispersion 
profiles are less accurately measured. 

There are problems and issues that cannot be addressed with a sample of this size. A much 
larger sample would likely find that a common value of [3 would not work for both red and blue 
galaxies, since it is expected that on the average blue galaxies will avoid the core, hence near the 
centre have smaller (3 than the red galaxies, and at large radius include many objects on their first 
cluster crossing, hence have higher (3 than the red galaxies. Furthermore the galaxies themselves 
are not invariant mass points. The blue galaxies are being altered by the cluster environment such 
that some of the their members are likely leaving the blue population to join the red population 
(Abraham et al. 1996). 
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The virial mass calculated from the full sample, empirically adjusted for its measurement 
biases, will correctly estimate the mass enclosed. With these two subsamples we have now 
tested each step in the chain of logic which supports our corrected, population adjusted value of 
fio = 0.19 ± 0.06. There is no compelling evidence for any remaining systematic errors of the 
cluster M/L as an estimator of the field value over the redshift range of this sample. 

We thank the Canadian Time Assignment Committee of the CFHT for allocations of observing 
time, and the CFHT organization for the technical support which made these observations feasible. 
Funding was provided by NSERC and NRC of Canada. 
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Fig. 1. — The projected number density profiles of the blue (solid squares and dotted lines) and red 
(circles and solid lines) cluster galaxies. The la confidence range from a Bootstrap error estimate 
is shown. 

Fig. 2. — The RMS line-of-sight velocity dispersion, (Tp j clS 3j function of projected radius. The blue 
galaxies (solid squares and dotted line) have a larger a p than the red ones (circles and solid line). 
The a p (R) are fitted with (3 = [0.0, 0.5]. 

Fig. 3. — The bias function refered to the total light profile, bMv{f) calculated from the mass 
profiles derived from blue (dotted line) and red (solid line) subsamples. The upper line for both 
subsamples is for (3 = 0, the lower line is for (3 = 1/2. 

Fig. 4. — The value of buv evaluated at r 2 oo- The blue subsample results are denoted with solid 
squares and the red with open circles. 
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